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TECHNICAL MEMORANDUM X-53268 


A METHOD OF IMPLEMENTING CUTOFF CONDITIONS FOR SATURN V LUNAR MISSIONS 
OUT OF EARTH PARKING ORBIT ASSUMING A CONTINUOUS GROUND LAUNCH WINDOW 


SUMMARY 


~ C i. tT O 


A method of implementing Saturn V lunar missions from an earth park- 
ing orbit is presented. The ground launch window is assumed continuous 
over a four and one-half hour period. The iterative guidance scheme 
combined with a set of auxiliary equations that define suitable S-IVB 
cutoff conditions, is the approach taken. The four inputs to the equa- 
tions that define cutoff conditions are represented as simple third- 
degree polynomials as a function of ignition time. 


Errors at lunar arrival caused by the separate and combined effects 
of the guidance equations, cutoff equations and input representations are 
shown. Vehicle performance variations and parking orbit injection errors 
are included as perturbations. 


Appendix I explains how aim vectors were computed for the cutoff 
equations. Appendix II presents all guidance equations and related imple- 
mentation procedures. Appendix III gives the derivation of the auxiliary 
cutoff equations. _ 


No error at lunar arrival was large enough to require a midcourse 
correction greater than one meter per second assuming a transfer time of 
three days and the midcourse correction occurs five hours after injection. 
Since this result is insignificant when compared to expected hardware 
errors, the implementation procedures presented are adequate to define 
cutoff conditions for Saturn V lunar missions. 


ODUCTION 


The iterative guidance scheme will generate steering functions which 
will insure attainment of a desired cutoff condition. Since these steer- 
ing functions are updated in flight, the cutoff condition does not have to 
be invariant with respect to time. Each evaluation of the guidance scheme 
generates steering functions which will direct the vehicle, in a near 
optimum maneuver, toward whatever end condition was defined for that 
evaluation . 


A set of equations describing a cutoff condition which will satisfy 
a given mission is defined as a hypersurface. The hypersurface used in 
this report requires four basic inputs: a unit aim vector, eccentricity, 
cutoff energy and an aim vector magnitude. For a continuous launch capa- 
bility, these four quantities are generated in the guidance computer as 
a function of parking orbit ignition time. The unit aim vector was imple- 
mented by fitting each component as a third-degree least square polynomial. 
The three magnitudes, eccentricity, cutoff energy and aim vector magnitude, 
were also represented as third-degree least square polynomials. 

All nominal trajectories used to compute basic input data for the 
hypersurface were optimized for cutoff weight at lunar injection subject 
to three lunar end conditions: (1) radius of closest approach equal to 

1923 kilometers (RCA) , (2) a flight time of 72 hours from lunar injection 
to radius of closest approach (Tf ) , and (3) inclination of the flight 
plane with respect to the moon’s equatorial plane at arrival minimized for 
the chosen launch time (INC) . For the particular launch day selected for 
the error analysis, the declination of the moon at arrival was near 
minimum. Therefore, launch azimuth is not symmetrical with respect to 
launch time. 

To evaluate the effectiveness of the method used to compute basic 
hypersurface inputs, implementation procedures, the guidance equations, 
and hypersurface accuracy, a set of error analyses is shown. A measure 
of the accuracy of any lunar injection guidance scheme is how far the 
nominal radius of closest approach (RCA) is missed. Each error analysis 
shows the RCA miss distance, time of flight variation, inclination errors 
and payload losses for six. launch azimuths ranging from 72 to 105 degrees. 
This corresponds to a launch window of about 41/2 hours. What each error 
analysis shows is briefly described as follows: 

(1) Errors and performance losses caused by the guidance 
scheme and hypersurface . 

(2) Errors and performance losses caused by the hypersurface. 

(3) Errors and performance losses caused by the guidance 
scheme, hypersurface and implementation of basic hyper- 
surface inputs by polynomial representation. 

(4) Same as (3) except vehicle performance and parking orbit 
perturbations are superimposed. Only a 72 degree launch 
azimuth case is shown, since it is representative of 
other azimuths. 

For reference, the calculus of variations optimum nominals are shown. 

From these trajectories the basic hypersurface inputs were calculated as 
described in Appendix I. 
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II. GUIDANCE SCHEME AND HYPERSURFACE ERRORS 


Errors at lunar arrival caused by the guidance scheme and hyper- 
surface equations are shown in Table I. Nominal input values for the 
hypersurface were taken from the optimized nominal trajectories shown 
in Table V. Initial conditions and vehicle characteristics were the 
same as the optimized nominals. 


TABLE I 


Azimuth 

(deg) 

ARCA 

(km) 

AT p 

(sec) 

AINC 

(deg) 

AW 

(lbs) 

72 

-3 

-4 

-.01 

-24 

77.7 

2 

9 

.00 

-25 

90 

-8 

i 

i 

-13 

-.01 

-25 

93.7 

i 

1 

-9 

-17 

-.05 

-25 

98.4 

i 

2 

8 

.00 

-25 

105 

-2 

-5 

.00 

-24 


The radius of closest approach miss distance (ARCA) is shown in 
kilometers. The deviation time of flight (ATp) from the 72-hour nominal 
value is shown in seconds. Inclination errors (AINC) are usually very 
small, and are shown in degrees. Payload losses (AM ) , due almost entirely 
to the guidance scheme, are shown as the number of pounds the guided runs 
were below the optimized nominals whose payloads were approximately 
128,000 lbs. at lunar injection. The same guidance equations which were 
used for ascent into an earth parking orbit were used for the lunar 
injection phase. Slight modification of the ascent guidance equations 
can reduce payload losses to about three pounds. However, since either 
case is acceptable, the same guidance equations are used for both opera- 
tions . 

Appendix II shows in detail the guidance equations, inputs required, 
and hypersurface equations used for these error analyses. 
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III. HYPERSURFACE ERRORS 


Errors at lunar arrival caused by the cutoff hypersurface equations 
only are shown in Table II. The powered flight phase was optimized for 
payload by calculus of variations techniques subject to the cutoff con- 
dition specified by the hypersur face . It is emphasized that Table II 
does not represent guided trajectories. 


TABLE II 


Azimuth 

(deg) 

ARCA 

(km) 

ATp 

(sec) 

AINC 

(deg) 

AW 

(lbs) 

72 

-3 

-3 

-.02 

0 . 

77.7 

-4 

-2 

.02 

0 . 

90 

-1 

8 

.16 

0 . 

93.7 

-3 

4 

.10 

0 . 

98.4 

-2 

-1 

.00 

0 . 

105 

-2 

-7 

• 21 

0 . 


Therefore, the deviations from the optimized nominals represent only hyper 
surface errors and cutoff tolerance errors. The cutoff tolerances were 
small enough to enforce RCA to within one kilometer and flight time within 
two seconds. Inclination enforcement was much more difficult to achieve, 
and the differences between the values shown in Table I and II show the 
guidance scheme can cut off in a given plane more accurately than a cal- 
culus of variations program using a reasonable number of isolation runs. 
Experience has shown that this is true for other types of missions as well 
Payload losses were less than one pound for all cases. 


4 













IV. REPRESENTATION ERRORS 


If the basic hypersurface inputs (M, M, e, c 3 ) are represented as 
third-degree polynomials as a function of ignition time (measured from 
midnight) , larger errors at lunar arrival will occur than if nominal 
hypersurface inputs are used. Table III represents guided trajectories 
with hypersurface inputs curve fit as a function of parking orbit 
ignition time. Ignition occurred when the vehicle was a fixed angle 
from the aim vector. Since the aim vector is a function of ignition 
time, a series of tests is necessary to determine ignition time. This 
will be no problem for the onboard computer, and the ingition criterion 
is satisfactory. 


TABLE III 


Azimuth 

(deg) 

AR.CA 

(km) 

at f 

(sec) 

AINC 

(deg) 

4 W 

(lbs) 

72 

-15 

9 

' .47 

-26 

77.7 

30 

45 

1.63 

-31 

90 

-45 

-83 

.29 

-25 

93.7 

-24 

-86 

1.64 

-30 

98.4 

8 

-51 

3.57 

-40 

105 

-24 

-56 

.46 

-24 


After time of ignition was determined, the polynomials were evaluated 
and the results used as inputs to the hypersurface. 

A comparison of Table III with Table I shows the radius of closest 
approach error increased slightly because of representation errors. The 
accuracy decrease is a small consideration when compared to the simplicity 
of implementation for an entire launch window. Larger errors are shown 
in this table than would be expected, because the launch window was con- 
sidered to be 4 1/2 hours and the maximum allowed for a realistic mission 
will be 2 1/2 hours. Therefore, the nonsymmetry of the hypersurface inputs 
versus time resulted in larger curve-fit errors than if the more linear 
part of the launch window were chosen. 
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V # VEHICLE AND PARKING ORBIT PERTURBATIONS 


Assuming the S-IVB vehicle will not perform nominally and that park- 
ing orbit injection conditions will not be perfect, Table IV presents 
lunar arrival errors caused by these type of perturbations for a 72-degree 
launch azimuth. The results are representative of any launch azimuth 
contained in the launch window. 


TABLE IV 


Perturbation 

• ARCA 
(km) 

AT f 

(sec) 

AINC 

(deg) 

A W 

(lbs) 

Nominal 

-15 

+9 

.47 

-26 

+F 

8 

-23 

.39 

-29 

-F 

-47 

-37 

.52 

-29 

+I sp 

-24 

-30 

.46 

-25 


-16 

-31 

.44 

-25 

+W 

-23 

-26 

.47 

-25 

-w 

-16 

-36. 

.42 

-25 

+R 

-3 

-36 

.41 

-25 

-R 

-38 

-25 

.49 

-25 

+V 

-16 

-25 

.45 

-24 

-V 

-11 

-9 

.49 

-27 

+t i 

-37 

-11 

.58 

-25 

-ti 

12 

-16 

.39 

-25 

+t L 

-34 

-70 

- .04 

-28 

-t L 

6 

32 

1.84 

-26 













A typical S-IVB vehicle was used for the nominal case. Magnitudes of the 
perturbations were as follows. 


±F 

±8000 

± ^sp 

±8.62 

±W 

±2500 

±R 

±30 

±V 

±15 

±fc i 

±15 

~ t L 

±1 


lbs thrust 

sec specific impulse 

lbs weight uncertainty 

km parking orbit altitude variation 

m/sec parking orbit velocity variation 

sec time of ignition error 

minute time of launch (azimuth mis -alignment 


error) . 


These magnitudes are far larger than any expected vehicle perturbations 
or parking orbit variations. 

If the cutoff ellipse were invariant even when a parking orbit alti- 
tude variation occurred, large payload losses would result. However, 
small payload losses are maintained by varying the eccentricity of the 
cutoff ellipse so that the altitude gain by the optimum nominal profiles 
is enforced. This is accomplished as a part of the cutoff hypersurface 
equations by computing the change in eccentricity as 


Ae 


_ /V' 

" w 


Ar 


where de/dr is an analytic first order approximation. The cutoff hyper- 
surface equations are shown in detail in the flow chart on page 16. 
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VI. OPTIMUM NOMINALS 


The problem of computing optimum nominal lunar transfer trajectories 
from an earth parking orbit requires a large amount of computer time. 

A calculus of variations powered program assuming constant thrust and 
flow rate was used. Once a parking orbit is established and the lunar 
end conditions specified, the problem can be resolved to finding the 
minimum time required to transfer the vehicle from some set of initial 
state variables to some other set of terminal state variables subject to 
the lunar end conditions. The initial state variables are only a func- 
tion of time of ignition. This problem can be solved in five to ten 
minutes on the IBM 7094 computer. 

Table V shows the nominal trajectories that were used to obtain 
input data for the cutoff hypersurface. Small variations in RCA and 
flight time are caused by isolation tolerances. Time of launch is shown 
to illustrate the nonsymmetry of launch azimuth versus launch time. A 
3-hour launch window is possible between 72 and 90 degrees launch azimuth. 
However, the hypersurface inputs were third-degree polynomials in time 
from 72 to 105 degrees launch azimuth. 


TABLE V 


Azimuth 

(deg) 

RCA 

(km) 

(hr/ sec) 

INC 

(deg) 

Time of Launch 
(hr/min) 

72 

1921 

72/7 

1.11 

6/25 

77.7 

1922 

72/5 

.32 

7/32 

90 

1922 

72/3 

.75 

9/43 

93.7 

1921 

72/3 

.70 

10/7 

98.4 

1923 

72/13 

.30 

10/27 

105 

1922 

72/3 

.46 

10/42 
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VII. CONCLUSIONS 


Lunar arrival errors caused by the iterative guidance scheme, cut- 
off hypersurface, and representation of inputs by simple polynomials are 
insignificant when compared to expected errors caused by vehicle hard- 
ware. Therefore, the implementation procedures as presented are adequate 
to define cutoff conditions for the S-IVB vehicle. 

No plane change cases have been presented. Current studies indicate 
that a slight modification of the hypersurface equations will increase 
performance for missions that require a plane change maneuver. 
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APPENDIX I 

Aim Vector Calculation 


After a nominal powered trajectory is obtained which satisfies all 
mission constraints, basic inputs required by the hypersurface can be 
computed. The basic inputs are a unit aim vector (M) , eccentricity (e), 
cutoff energy (c 3 ) , and an aim vector magnitude (M) . An optimum nominal 
trajectory is not essential. However, the set of equations that define 
the hypersurface will assume the cutoff condition of the nominal only 
when furnished the nominal range angle. Since the range angle as com- 
puted by the Iterative Guidance Scheme* will be nearly optimum, if the 
reference is not an optimum trajectory, nominal cutoff will not occur. 
Therefore, it is usually desirable that the reference trajectory be an 
optimum so that any deviation in cutoff conditions will cause performance 
losses . 

All trajectories used to define basic hypersurface input values 
were optimized for payload subject to mission constraints. No plane 
change maneuver was performed. Payload was optimized subject to an 
energy cutoff value that resulted in a desired flight time. Time of 
coast and time of launch were varied to meet the other constraints 
described in the introduction. The maximum payload possible for a given 
flight time was computed. 

Using these optimum nominal powered flight trajectories, the basic 
hypersurface inputs are computed as follows with all vectors computed 
in the space-fixed plumbline’ earth-centered coordinate system. 


r = 
v = 
r = 
v = 

N = 


position vector at cutoff 
velocity vector at cutoff 
magnitude of r 
magnitude of v 
r x v 


r x v 


- 


7 2 _ lli ^ = .3986032 x 10 15 m 3 /sec 2 
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e = (1 + c 3 r 2 v 2 sin 2 0/ja 2 ) 1 ^ 2 


a 


SJl 

M* 


r = N x r 

~p - - 


0 = arc cos < — 

c I e 


- (1 - e 2 ) 
r y 


After powered cutoff the vehicle coasts to the desired lunar end 
condition. Any point along the ballistic flight could be chosen as the 
aim point. For all error analyses in this report the aim point was 
chosen where the vehicle entered the lunar sphere of influence (defined 
as 35,000 statute miles from the center of the moon). The distance from 
the center of the earth when the vehicle entered the lunar sphere of 
influence was chosen as the aim vector magnitude. 


M = aim vector magnitude 


0 = arc cos < — 

T ] e 


5 (1 - **> - 1 


M = - 
— r 


r cos (jZL - 0 ) + r sin 
— V/ T 'c' — p 


Now the four basic inputs required for 
computed: 


} 

<0 X - 0 C ) ■ 

the hypersurface have been 


M = unit aim vector 
M = aim vector magnitude 
e = eccentricity 
c 3 = cutoff energy. 
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It is important to select the aim point from an integrated ballistic 
flight that satisfies a desired end condition. The desired end condition 
need not correspond to lunar arrival conditions even though that is the 
final objective. For example, an aim point may be chosen where a mid- 
course maneuver is to occur. 

Reviewing the equations used to compute the hypersurface inputs 
will reveal why lunar arrival errors are small when cutoff is approxi- 
mated by simple elliptic equations. The desired cutoff energy and 
eccentricity are computed using the cutoff condition of the nominal. 

An arbitrary distance from the center of the earth is chosen as the aim 
vector magnitude. This magnitude determines the angle between the aim 
vector and perigee vector. Next the aim vector is computed in the 
nominal cutoff plane. Although the actual cutoff plane as determined 
by the hypersurface will depend upon parking orbit variations, the hyper- 
surface equations will insure that the injection point defined by the 
hypersurface will be near the injection point of the nominal. This means 
that both ballistic trajectories will experience almost the same gravi- 
tational perturbations. The small differences will cause only small lunar 
arrival errors . 
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APPENDIX II 

Guidance Equations and Hypersurface 
A. Implementation Techniques 

The guidance equations and hypersurface used for the error analysis 
shown in the text are presented in the form of a flow chart on page 16. 

Ignition occurred when the vehicle was a fixed angle from the aim 
vector (TIGA was set equal to zero). This ignition criterion was used 
for all cases where the hypersurface inputs were represented as poly- 
nomials. 

The hypersurface equations themselves can be considered the equa- 
tions necessary to compute the perigee vector (S) . All other equations 
are used only to translate from the hypersurface coordinate system to 
the guidance reference system. 

The guidance equations themselves begin in the block where the 
terminal range angle is estimated and continue until the pitch and yaw 
angles are computed. 

In addition to the equations used to compute pitch and yaw, the 
computer simulation program used a pitch and yaw rate given as 



and 

K 3 A 

"T 1 * 

y 

respectively. 

Therefore, the equations for Xp and Xy (chi pitch and chi yaw in the 
guidance reference coordinate system) become 

X' = X - K, + K 2 St 
P P 1 

Xy = Xy - K 3 + K 4 St, 

where small St is defined as the time since the guidance equations were 
last evaluated. 
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ITERATIVE GUIDANCE MODE EQUATIONS FLIGHT OUT OF ORBIT 


INPUTS FOR OUT OF ORBIT 


Vehicle Dependent Inputs (!) 

T* * estknoted burning time out of 
parking orbit 

Mission Dependent Inputs (28) 

M. »M |# + M n T+M |t T*^M, s T s 

unit wn vector 
M •R 0 +R,T^R # T 2 + R s T* 
aim vector magnitude 
e H -e<)+ei T+e t T 2 +e 3 T 3 
nominal eccentricity 
C 3 *C 0 + C| T + C Z T*+C 3 T 3 
cutoff enery 

r„ « nominal radius ot time of ignition 
K, ■mission constant 
COSA • cosine of angle before M to ignite 
TiGA ■ cosine of angle before M to begin 
chill down 


FiFl-lsl 


t 

AFTER PARKING ORBIT INJECTION 
COAST UNTIL 

COS A « ♦TIGA 

IGNITION OCCURS IN TIG SECONDS 


MATRIX TO TRANSFORM FROM PLATFORM 
COORDINATES TO GUIDANCE COORDINATES 


w* 1 

sin A z 

0 

cos A z 

0 

1 

0 

-cos A z 

0 

sin A z 


1 

0 

0 

[■]- 

O 

cos* 0 

sin 

L ° 

-sinfc, 

cot to 

[Cl- 

cos e 

-sine 

0 

sine 

cose 

0 

0 

0 

l 

io]. 

cosi 

0 

sini 

0 

1 

0 

-sini 

o 

cosi 

[•]-[»] [c] 

[•1 [*] 



COMPUTE AT SECOND S-IVB IGNITION 


M unit aim vector 
M aim vector magnitude 
e N nominal eccentricity 
C 3 cutoff energy 
e 4„<e N -l)*l 


eccentricity of cutoff ellipee 
semilatus rectum 
true anomaly of aim vector 
constant 

normal to parking orbit plane 


K 5 »>/u7p 

DD«(N constant 

FF«N M/INXMI* constant 
§ « [dD M- FF W + - DD)i] cost* perigee vector 

S A ■ vector perpendicular to perigee vector 

Qi *§x§ A normol to cutoff pkme 

fl -cos A 2 cost,, l+sint 0 sin A z cos ^ k earth spin vector 
cos i * g Q, inclination of cutoff plane 
I L *-ooe A z sin fl,J+cos£ 0 j 

♦ ski A z sin A 0 fc vector to launch site 
ON * £i x & constant vector 

8 ■arc ton descending node 
' I I ' 

a *arc tan ( § f j f 1 ) frofln § t0 <*®*c*nding node 


TERMINAL RANGE ANGLE CALCULATION 
IN ORBIT PLANE 

L 2 *ln[r 2 /(r 2 -T 2 )] 

*1* * v en 2 ^ T 2L2“T 2 ) 

S e * J*~T z v ex z L-2 
Q*-iVe Xl T 2 l + r 2 S 2 
t, - V„,/(F/M) 

L| * ln[r, /(t,-T,)] 

Ji * V M( (r,L r T, ) 

Si ■ Jr T, V #B| Li 

Q. = i Vm, Ti 2 +t, s, 

L»* V M( L| ♦ V„ t L* 

8i ■ -Si - Jt ♦ L*( T 2 ♦ T* ) 

Tic* T|+T 2 +T c 
L' » ln[v(r 5 -T')] 

J* * V„, <t 3 L'-T*) 

T*- T'+ T, c 
A' ■ L* + V„ S L' 

8*» VTh-j'+A'T'-K.Ua-T’HA'+V-V^L* 


*fp 

■M 

V 


y 

*ff 


z 


tr* arctan(x FP /yp P >+ ^ (8, + ®*) 

Ad> ■ t? + cc true anomaly 

«T»«rc..n[ | ; e ,l c " 0 ^ ] P0lh. n ,l. 

^ T * p/(l + eco*A4>) radius 

V T " K 5 (l+e 2 + 2ecos Ad>)* velocity 

V T sin 6j velocity along if axis 

( T « v T cos 8 t velocity along ( axis 


POSITION, VELOCITY, F/m a 
p. FROM NAVIGATION 

x, y, i, x t y, i, F/m, ft 


COORDINATE ROTATIONS 


Ysr* 

M* 
M • 


"cotfr -slntr 
sin^ T co*t T 
<5 o 

MM 


si 



x 

& 

X 

■W 

y 

* 

7 • [k] y 


z 

u 

ij L f . 



s 


TIME TO CUTOFF 

CALCULATION 

as*, -t-i, T* 
,r(A<«)* + (AV 

)*+(acv .n 

6 -*L A' 

AT *G(t b -T')/V„ 3 

t 3 -t'+at 

T* ■ T*+ AT 

-Aj 


1 


PITCH AND YAW ANGLE CALCULATIONS 

, 

M ■ jcos x y + K 5 sinX y 

A i’ Ai*-(S,AT 

N ■ 1 K4 sin J(y 

An* a^-b^at 

Ap ■ MAy - NBy 

At* At*"tjAT 

p l ■-|T| 2 V0K,+ T,J| 

X y 9 ore ton[AC/(A^ 2 +A^ 2 ) ,/2 ] 

F*-iT* 2 Vi, 2 + 4 i(Ti+ 2 T|» 

5( p * arc tan ( An /AC) 

P 3 * * i T 3 2 Vsxs + J 3 < r 3 + 2 T| C ) 

VsxjLj * W‘ +G 

B P * MB, - N(T| 2 V #K2 L 2 + 

A,* L«*V„ S L, 

+ T,|vi,,Lj + 

J) ■ J'+ TjG 

+ 1} + P 2 + Pj) 

C* ■ J3 - AyTj - 8| 

e p* n~nT + nT*+ ingT* 2- 

Cy ■ c‘ COS Xy 

-(MC'-ND')tin 

£,. C+{T*+iC,T’ 2 -C'.irx, 

Cp- (MC'-ND')cosxp 

83* J| + J2 + T,V #X2 L 2 

o'. |T,*V„ s *t,S s 

By * J 3 + T IC v ex 3 l- 3 + 83 

«l* iVVo.+ r.O, 

84* Q| ♦ O2- T2J1+ TjS2 

Us- iTj J V„j+< 32 (rj-<- 2 T|) 

8 5 * 83 (T s ♦ T c ) 

U J * iVv„ J + 0 , (rj + 2T le ) 

S 3 * J 3 “ T 3 V tX 3 L 3 

Dp* {mo'-n[u, + u 2 + u s + t, j S2 + 

oj- ±ifv,*WV 

♦T|cSj-(Tj+T c )(P,+P 2 + 

D* ■ 84-83+03 

+ T l 2 V„ i L2)-P 1 T2]}cos Xp 

Dy- D' COS Xy 

K| ■ Bp Ep/ ( ApDp - Bp Cp) 

K3 ■ ByEy/(AyDy-ByCy) 

Xp* Xp- K . 

Xy* Xy" K 3 


CUTOFF WHEN 



















It is important to the guidance implementation procedure to fix 
the end conditions as predicted in the terminal range angle block about 
25 seconds before cutoff. Normally K x , K 2 , K 3 and K 4 are frozen 10 sec- 
onds before cutoff. This procedure prevents large turning rates from 
building up after the guidance parameters are frozen. 


To rotate the direction cosines of the guidance coordinate system 
to an earth-centered plumbline coordinate system, the inverse of the K 
matrix is used. This is the system required in the flight computer. 

Note that although the guidance scheme computes* chi yaw before 
chi pitch, the rotation order is chi pitch then chi yaw, which is the 
same as the vehicle platform. 


— 


— — 

DCX 


cos X cos X 

p y 

DCY 

= 

sin X cos X 

p y 

DCZ 


- sin X y 

_ T _ 
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B. A Description of the Iterative Guidance Scheme 

Some of the basic principles of the iterative scheme can best 
be demonstrated by assuming constant gravity, thrust and specific impulse. 
If only a velocity end condition is enforced, it can be shown that a 
constant thrust direction is the optimum steering law [1], The thrust 
direction can be found by the following geometric construction. 



I 
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FIGURE 1. THRUST DIRECTION. GEOMETRY 






If the remaining time, T*, is not known, it can be estimated by 
the characteristic velocity equation and updated for each guidance evalua- 
tion cycle. Terminal velocity components and yx are assumed constants. 
As time-to-burn approaches zero (cutoff), the steering commands are frozen 
to avoid indeterminate expressions. Since actual vehicle implementation 
procedures require frozen steering functions near cutoff, no scheme 
generality is sacrificed. 

A constant thrust direction cannot enforce both altitude and 
velocity. Therefore, a different steering law is needed. The steering 
law which the iterative uses is 


X = X - K x + K 2 5t, 


where K x and^K 2 are subject to the conditions that the velocity constraint 
enforced by X is not violated. Each evaluation cycle updates X, and K 2 
using the current state variables and vehicle characteristics. Small &t 
is the time since the last guidance evaluation cycle. 

This new steering law has a theoretical basis: the first order 

expansion of the calculus of variations steering law derived subject to 
the same assumptions (flat earth, constant thrust and specific impulse, 
altitude and velocity enforced, no range enforcement) is 


X = a + bt 


where a and b are computed subject to current state variables, vehicle 
characteristics and desired terminal conditions. 

It is possible to avoid specifying lateral terminal displacement 
and velocity components. This is accomplished by defining a guidance 
reference coordinate system that has one axis perpendicular to the final 
cutoff plane (Figure 2). All calculations by the guidance equations are 
done in the guidance reference coordinate system. The coordinates are 
defined as follows: r\ is measured from the center of the earth and 

passes through the final cutoff point, | is measured from the center of 
the earth, perpendicular to rj, in the direction of flight, and £ completes 
a right-handed coordinate system. 

Usually, cutoff conditions are specified as a desired radius (r), 
velocity magnitude (v), and path angle (e), measured from the local hori- 
zonal. These cutoff conditions transformed to the g, q, £ system are 
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FIGURE 2. GEOMETRY OF GUIDANCE REFERENCE SYSTEM 
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i T ” 0 


i T = v cos e 
T} T = r r| T = V sin 0 

C T - o E T - o. 

To rotate platform coordinates to the desired cutoff plane, 
the G matrix is computed by the launch -vehicle computer. This matrix is 
a constant and can be computed as soon as launch azimuth, launch latitude, 
descending node and inclination of the cutoff plane are known. The 
terminal range angle, is the remaining rotation necessary to transform 

from the cutoff plane to the guidance reference system. The matrix that 
rotates the platform coordinates into the guidance reference system is 
called the K matrix. It is computed in the coordinate rotation section. 
The K matrix is updated through every evaluation cycle so that rj 
always passes through the predicted cutoff point. Notice that the pre- 
dicted terminal values are in the guidance, reference system so that 
and are zero (| T is zero also). 

Time-to-burn is updated based on velocities to be gained measured 
in the guidance reference system. Cutoff occurs when the vehicle reaches 
the desired energy. At this time the altitude and path angle constraints 
will be enforced. 
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APPENDIX III 


Derivation of the Cutoff Hypersurface Perigee Equation 


It has been empirically shown [14] that the injection points of 
earth-moon trajectories with constant transit times map concentric 
circles on the surface of the earth. The diameter of these injection 
rings is a function of the injection path angle. All nominal optimized 
trajectories generated for a particular vehicle from a circular parking 
orbit will have almost equal injection path angles. The major lack of 
path angle uniformity is caused by earth gravitational variations result- 
ing from different launch azimuths. Therefore, the injection ring can 
be considered to have a constant diameter for the type of trajectories 
considered in this report. 

Each injection point defines a transfer ellipse. The perigee points 
of these transfer ellipses map a perigee circle which is smaller and con- 
centric to the injection circle. For convenience of representation, the 
cutoff hypersurface equations are derived based on the perigee ring. 

Consider one particular optimum transfer trajectory where an aim 
point has been selected and the unit aim vector computed as described 
in Appendix I. A rotation of the trajectory plane about the aim vector 
will cause the perigee vector corresponding to the transfer ellipse to 
sweep out a right circular cone. The apex of this cone is the center 
of the earth and the perimeter of the base is the perigee ring. Notice 
also that the axis coincides with the unit aim vector. 

Since all optimum nominal cutoff planes contain the axis of the 
perigee cone, the cutoff plane for all earth-moon trajectories is defined 
to be the plane that contains a perigee vector in the parking orbit plane 
and the aim vector. The perigee vector in the parking orbit plane will 
be the vector from the center of the earth to the point of intersection 
of the parking orbit plane and the perigee circle ahead of the ignition 
point. Actually, both the perigee and aim vectors are considered unit 
vectors with separate magnitudes. 
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M unit aim vector 

unit perigee vector 
P angle between M and S_ 

N unit vector normal to parking orbit plane at ignition. 


For clarity of derivation, the aim vector is assumed to point to 
the center of the perigee circle, and the acute angle between M and S_ 
is defined as p. The perigee vector is computed in an orthogonal 
coordinate system defined in terms of the two unit vectors N and M. 
The axis of the coordinate system are 
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M 


unit aim vector 


N x M 
' N x M ’ 


unit vector in the parking orbit plane 


M x (N x M) 
_ Nx M 


unit vector completing the right-handed coordinate 
system 


where N and M are computed in the earth centered plumbline coordinate 
system. Since any vector can be expressed as a linear combination of any 
other three non-coplaner vectors, let 


S = a M + b 


M x (N x M) 
N x M 


+ c 


N x M 
' Nx M 


where a, b and c are the direction cosines of £ to be determined. By 
observation 


and 


a = cos p 


b = - cos p tan 7, 

where 7 is defined by 

sin 7 = N • M 
cos 7 = |N x M| . 

Now c can be determined by 

c = \ll - a 2 - b 2 = \l 1 - cos 2 p sec 2 7, 
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or, in terms of p, N and M, 


c 


f\ 


N x M 


COS 2 p 


N x M 


Since 


M x (N x M) = N - (M • N) M, 
S can be written as 


S 


b 


M • N \ N 

t— rrr I M + b -p— ; rr 

N x M y - NxM 


N x M 
N x M| * 


Substituting a, b, c, -M for M and p = jt - the equation for S_ becomes 


S 


DD M 


FF N + 


NxM 
N~ x M 



cos 0* 


where 


DD 


(N • M) 2 
‘n x M| 2 + 1 


FF 


N • M 
N x M I 2 


This is the form of the perigee vector used in the precompute section on 
page 16. 
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A unit vector normal to the cutoff plane can be written as 


S x M 
s i n 


Using the normal to the desired cutoff plane, the inclination and descend- 
ing node can be computed. Next the location of the perigee vector with 
respect to the guidance reference coordinate system is computed. The 
cutoff path angle, radius and velocity are computed as functions of true 
anomaly of cutoff, which is updated every guidance equation cycle by the 
guidance scheme. 
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